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Abstract. We propose a scheme to simulate the dynamics of neutrino 
oscillations using trapped ions. For neutrinos in 1+1 dimensions, our scheme 
is experimentally implementable with existing trapped ion technology. We show 
that the three generation neutrino oscillations can be realized with three ions for 
1+3 and 1+1 dimensions where the latter case only requires experimentally proven 
two-ion interactions. For this case, we discuss two setups utilizing different types 
of spin-spin interactions. Our method can be readily applied to two generation 
neutrino oscillations requiring fewer ions and lasers. We give a brief outline of a 
possible experimental scenario. 
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1. Introduction 

Ever since Pauli inferred their existence, the exact nature of neutrinos have been 
a mystery due to their tendency to avoid interacting with other particles. One 
prominent question was whether the neutrinos are massive particles. In 1957 
Pontecorvo suggested that massive neutrinos can oscillate (change its flavour) and 
later noticed that it could explain the solar neutrino problem [1]. Since then large 
scale experiments have confirmed that neutrinos have small but nonzero masses and 
give rise to oscillations between fiavour eigenstates [5]. The discovery of neutrino 
oscillations is regarded as one of the most important discoveries of modern elementary 
particle physics and have lead to many proposals on physics beyond the standard 
model. 

Recently there has been a growing interest in simulating exotic relativistic 
phenomena using other controllable quantum systems. There have been proposals to 
simulate equation in curved spacetime [3], Unruh effect [4], and black- hole properties 
in BECs [5] . More recently, a scheme to simulate the many-body dynamics of a Dirac 
particle, viz., Schwinger effect, has also been proposed [6]. The seminal work on 
quantum simulation of the Dirac equation using trapped ions by Lamata et al. [7] and 
its subsequent experimental realization by Gerritsma et al. [5] are of particular interest 
to the current work. More recently, there has been a proposal [9] and experimental 
demonstration |10| to simulate the Klein paradox. 

On the other hand, trapped ion system is also one of the leading candidates for 
simulation of quantum spin systems. The idea was originally proposed by Porras 
and Cirac [llj and recently similar schemes were verified experimentally. These 
experiments include the simulation of a quantum Ising model with two |12j or more |13| 
ions, where phonon- mediated spin-spin interactions are realized [T4l[T5lfT6l[T7l[T8lll9) . 
Combining the interactions giving rise to the linear momentum term used in the 
simulation of the Dirac equation and the spin-spin interactions used in the simulation 
of the Ising models simultaneously broadens the types of systems that can be simulated 
with trapped ions. For example, Casanova et al. have considered quantum simulation 
of the Majorana equation and unphysical operations |20| . 

In this work, we show that neutrino oscillations can be simulated in trapped 
ion systems. Especially, observation of neutrino oscillations in 1+1 dimensions 
requires only the trapped ion technology demonstrated in experiments. Because 
whether the neutrinos are Dirac or Majorana particles bear no observable consequences 
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in neutrino oscillations, wc assume, without loss of generality, that neutrinos are 
Dirac particles as described by the minimally extended standard model. In this 
model, the charged current interaction creates a flavour eigenstate which is in a 
superposition of mass eigenstates. Our scheme paves the way for an experimental 
study of neutrino oscillations with controllable creation of the initial state and 
oscillation length, allowing different types of neutrino oscillations experiments (see, 
for example, |22j about different types of neutrino oscillations experiments) in a single 
setup. Also, non-trivial initial states not observed in nature can be readily created. 
For example, a state in a superposition of positive and negative energy eigenstates 
produces Zitterbewegung-like high frequency oscillations [23]. 



2. Standard theory of neutrino oscillations 

Here, we reproduce the gist of the standard theory of neutrino oscillations where 
the neutrinos are assumed to be created in a momentum eigenstate |22] . So far, 
experiments have verified that there are three flavours of neutrinos: electron, muon, 
and tauon neutrinos, which we denote as Uq. with a — e,fi,T. These flavor states 
are not mass eigenstates and therefore do not follow the dynamics given by the 
Dirac equation. They, however, are related to the mass eigenstates Vk by a mixing 
matrix through the equation \i/a) = J2k^ak\^k), where k = 1,2,3 is used to label 
different mass eigenstates. Since a massive neutrino state obeys equation, one can 
write {vkit)) = e~''^'=*|i^fe), where Ek = ±^y{c\p\)^ + {rrikC^Y, given that the state 
is in the momentum eigenstate with momentum p. The time evolution of a flavor 
eigenstate is given by 

k 

= T.T.u*^^^'''''''uM^ (1) 

P k 

which means that the probability for the flavour of the neutrino to change from a to 
/3 is 

In the ultrarelativistic limit c|p| 3> mc^ 

Ek-E,^ (3 

^ 2E ^ ' 

with Am^j = to| — mj and E = c|p|. Thus, after replacing the propagation time t 
with the distance traveled L/c, the probability becomes 

P..^., (t) = E K^UpkU^,U*^^e-^^^^. (4) 

Measuring the probability for flavour changes thus allows one to gain information 
about the squared mass difference and the mixing matrix. Especially, if all the masses 
are equal there would be no neutrino oscillations 
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3. Trapped ion implementation 

Here, we show how the three generation neutrino oscillations in 1+1 D can be 
implemented in a system of 3 ions utilizing only previously tested ion manipulations. 
A similar implementation, extending the original scheme in [7] for 3+1 D is also 
possible, but as this scheme requires a slightly more complicated experimental setup 
not realized in the lab thus far, we focus on the 1+1 D case. The 3+1 D case is briefly 
discussed later. 

3.1. 1+1 dimensions 

In 1+1 D the Dirac equation reads [M] 

= {cpiyx + mc^(T^) il>, (5) 

which can be simulated by a single trapped ion with two internal levels [7j [8] . The 
ion is driven by a bichromatic laser that couples the internal levels with motional 
sidebands to create the Dirac Hamiltonian 

Hd = 2rjAn<7xp + fla^, (6) 

where A = ^Xjlmuj is the size of the ground state wave function with ion mass m; uj is 
the frequency of the vibrational mode coupled to the internal states via the bichromatic 
laser. 77 is the Lamb-Dicke parameter, p is the momentum operator for the phonon 
mode, and the f2 term arises from the detuning 2i7 between the bichromatic light field 
and the qubit transition. This Hamiltonian is equivalent to the Dirac Hamiltonian in 
1+1 D with c = 2r/fiA and mc^ = VL. 

Neutrino oscillations arise from an interference between different energy 
eigenstates. To mimic this effect we need the correct relativistic energy dependence 
and an ability to create a superposition of different energy eigenstates, i.e. for ^ to be 
in a superposition of different mass eigenstates. For 1+1 D, each neutrino is described 
by a two-component spinor related to positive and negative energy states, which means 
that 6 basis states are required for 3 generations. We construct these basis states as 
follows: 

I, ^ _ (A999) \ \„x f "l55e) \ \„x ( "Ie53> \ (y. 
1"^^ - 1, P\geg) ) ' '^^^ " (, P\gee) ) ' l^^) " p\eeg) ) ' 

where \g) and |e) denote the two internal states of a qubit. It is easily seen that 
1 (g) ax 1 \vk) = Wk), where the ax on the r.h.s. exchanges the two basis states 
that define ji^fc). Also, it is easy to work out that 

Hss Wi) = {^1+^2) a, \i^i) , 

H.sWs) = {-ni + n2)a,\i^3), (8) 

where Hgs — fliaz <E} az <E) 1 + ^2^ 'i^ cr^ <E) az- Therefore with Eq. © we see that the 
Hamiltonian 

H = 1<^Hd<E) l + Hss, 

= 277Afi (1 (K) CTa, (X) l)p - ill (g) (Tz 1 

+ niaz(E)az<E)l + ^2i<S)az(S)az. (9) 
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produces the correct dynamics for three generations of neutrinos. That is, the mass 
eigenstatcs follow the Hamiltonian 

H\vk) = {caxP + mkc^az) \iyk), (10) 

with neutrino masses mic^ — n + + ^2, m2(? = + fii — ri2, msc^ = Q, — Q,i + Q,2- 
As mentioned earlier, the Hjj term can be created by focusing a detuned 
bichromatic laser on the second ion. The remaining part requires two-qubit gate type 
interactions on the qubits (1,2) and (2,3), respectively. One can get each of these two 
by selectively shining a pair of ions with a laser tuned to a particular phonon mode 
to mediate the ion- ion interaction [TH [151 UHl HZ] • The two sets of lasers should act 
on two different normal modes to avoid interfering with each other. For example, the 
two gate type interactions can utilize the center of mass mode and the zigzag mode 
in a transverse direction, while the linear momentum part utilizes the axial mode. 
Then, for the Hamiltonian to work for all times, one needs r/ <C 1 to avoid exciting 
a significant number of phonons which would result in qubit-phonon entanglement. 
Otherwise one could consider this as a gate operation which works only at certain 
times, in which case it is important to make sure that the two two-qubit interaction 
terms have commensurate gate times. 

Recently, Kim et al. have managed to create a tunable spin-spin couplings between 
trapped ions without excitation of real phonons [131 [T7] . This scheme offers nearly 
ideal spin-spin interactions that can be used to simulate neutrino oscillations. In their 
scheme all the ions are addressed simultaneously with two bichromatic laser beams 
whose optical beatnote detuning is far from each normal mode compared to that 
mode's sideband Rabi frequency. Thus the phonons are only virtually excited and the 
qubit states evolve according to the Hamiltonian 

at all times. Furthermore, the signs and magnitudes of Ji and J2 can be controlled 
by changing the beatnote detuning and the spin-flip Rabi frequencies, which in turn 
allows one to control the effective neutrino masses. To make use of this Hamiltonian, 
the basis states Iff), |e) should be changed to the x-basis states which we denote as 
|0), |1), and ions 1 and 2 should be swapped, i.e. 

/ a|000) \ / a|001) \ / a|010) \ 

- [ /3|100) ) ' l''^) - ^1101) ) ' I'^s) - ^1110) ) ■ (12) 

To mimic the dynamics of a neutrino, we add the linear momentum term and three 
extra single qubit lasers: 

H = 277Af^ {l®ay®l)p 

+ Jl {cTx ®(Tx®l + l®Crx® <7x) + J2<^x ®1®(Jx 

+ Jil®(Tx®l + Jil®l®ax- J(Jx 10 1, (13) 

yielding 

H\vi) = (2i]An(Typ + (7z{J + Ji + J2)) Wi) - JiWi), 
H\v2) = (2riAnayP + <Jz {J + J, - J2)) ^2) - J1W2), 

Hli^s) = {2r,A(layp + a, ( J - Ji + J2)) l^/g) - Ji\v^)- (14) 

Ignoring the constant term Ji, we get the Hamiltonian that describes 3 types of 
neutrinos with toic^ = J + Ji -V J2, m2C^ — J + Ji — J2, and m^c^ = J — Ji + J2- 
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Note that for the momentum term we now use an alternative but equivalent form cr^p, 
which can be implemented in the same way as the original scheme by changing the 
phase of the laser; also, we have assumed that the detuning is zero 




Figure 1. A potential experimental setup to simulate neutrino oscillations. The 
lasers Hi and Q,2 shine on all 3 ions with a given beatnote frequency and a 
wavevector difference Afc in a transverse direction, to create the mass terms v^ith 
the help of an extra laser acting on each ion (not shown on the figure). acting 
on the second ion, creates the linear momentum term of the Dirac Hamiltonian. 

The experimental setups of Islam et al. [13j and Kim et al. |17| use transverse 
phonon-modes to mediate the spin-spin interactions which means that the momentum 
part could utilize an axial mode or the perpendicular transverse mode. Then, three 
extra lasers with appropriate intensities and phases help create the appropriate mass 
terms. Figure 1 shows a possible experimental setup to simulate neutrino oscillations; 
for visibility three extra lasers addressing each ion are not shown on the figure. 

We also note that the rotated version of Eq. (|9]) (so that the gate operations 
Uz ® (Tz are replaced by Uy ® <jy for example) can be directly implemented using 
spatially dependent Rabi frequencies as opposed to a uniform Rabi frequency used in 
deriving the Hamiltonian in Eq. ([TT|) [2Tj . 

For two generations, only two trapped ions are required, which is interesting not 
only because they have less parameters, but also because some experiments are not 
sensitive to all three generations and can be described by an effective model with two- 
neutrino mixing j22j . The two generation case is significantly simpler than the three 
generation case as only one spin-spin coupling term and one single-qubit rotation term 
are needed to generate the mass terms. 

So far we have proposed schemes that utilize multiple trapped qubits. However, 
if one could find a stable multi-level trapped ion system, it is possible to simulate 
neutrino oscillations. For example, if there are 3 ground states and 3 excited states 
and a single bichromatic laser addressing the three transitions, one has the Dirac 
Hamiltonian ^ for each transition. If the transitions are not of the same frequency 
the masses would be different by default, otherwise one could use external fields to shift 
the energy levels. Possible difficulties in such single-ion schemes are short decoherence 
time and preparation of a general initial state. 
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3.2. 3+1 dimensions 

Here, we give a brief description of liow to simulate the full 3+1 D dynamics 
using the 4-lv scheme introduced in [7]. We can get the Dirac Hamiltonian for 3 
generations by considering the following basis states: = (|a), \b), |c), \d))®\a)®\a) , 
IV'2) = (|a), |&), |c), \d)) ® \h) ® \a), jV-s) = (|a), \b), |c), \d)) <g> \a) (g) \b). Apart from the 
Hamiltonian proposed in [7] acting on the first ion, we add the following spin-spin 
interaction terms -1^i(ct^'= - cr^'^) » «) 1 - i^2{cr^'' - cr^y'^) (S)a''/(g)l 
— r23(cr,y^ — ay"^) (E) 1 (E) cr^'' to obtain the correct rclativistic equations with mic^ = 
(f2 + ill), m2C^ = (ri + ^2), and m^c^ ~ (O + Jli + f^a), as can be easily verified. In 
most of the realistic cases there is only a single non-zero momentum component which 
reduces the number of lasers needed. 

4. A possible experimental scenario 

Once the Hamiltonian is engineered an appropriate initial state has to be prepared to 
simulate neutrino oscillations. For example, if the initial state is in a mass cigenstate 
one would observe no oscillations. These cases are however quite special and flavours 
would oscillate for a generic initial state. A physically interesting case is when the 
initial state is in a definite flavour state and has a momentum wave packet with a 
narrow momentum distribution around an average momentum in the ultrarelativistic 
regime. 

We describe a scenario where an electron neutrino is created which propagates 
for a certain amount of time before it is detected. For concreteness we use the 
implementation that uses cr^ ® iz type interactions. An electron neutrino state can 
be written as 

^ ^(V2\Mp)) ~W2{P))) , (15) 

according to a mixing matrix called tribimaximal mixing matrix that is consistent 
with experiments [25l [J^ . Note that this state only involves 2 states and thus can 
be described by a two generation model, i.e. with two ions. Since this state is not 
entangled: 

K{P)) - Iff) ® («p|.g) + Pp\e)) ® i (V2\g) - |e)) , (16) 

it can be prepared by performing single-qubit rotations on the ground state \ggg). 
Other flavor states are entangled and need 2 qubit gate operations. However, it has 
been shown that arbitrary states can be created by repeated use of gate operations 
[28] . and one such algorithm for trapped ions has already been proposed [l^. A 
similar proposal to prepare hadronic states in terms of up and down quarks and 
their spins has been given recently [SO]- The next step is to engineer the state of 
the phonon mode (momentum distribution). To create a state with a given average 
momentum, one would cool down the ions to the ground state and apply a momentum- 
displacement operation [8|. An arbitrary state can in principle be constructed 
by performing state-dependent displacement operations; positive or negative energy 
eigenstates have asymmetric spinor components that depend on the momentum of 
the particle and a superposition of momentum eigenstates gives a physical wavepacket 
of the particle. However, it is possible to approximate an energy cigenstate with a 
Gaussian momentum distribution for a Dirac particle in 1-1-1 D by approximating the 
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eigenstatc with a spinor that has an asymmetry between the average momenta of the 
components [S]- These approximate states can be created by focusing a momentum 
displacement laser on the second ion. Note that the asymmetry in the momentum 
distribution between the spinor components becomes smaller as one increases the 
average momentum of the particle, so in the ultrarelativistic regime the momentum 
wave function approaches the symmetric spinor ij^iji) oc exp[— (p — po)^/2(t](1, 1) with 
average momentum po- After creating the initial conditions the interactions can be 




2 4 6 8 10 12 14 
t(ms) 



Figure 2. Neutrino oscillations of an electron neutrino in a momentum 
eigenstatc. The top (black) curve represents electron component whereas the 
bottom (blue and red, overlapping) curves represent muon and tauon neutrino 
components. The kinetic energy is 27rx40 Khz and the rest mass energies are 
27rx(5, 6, and 7) Khz. 

switched on and the required states would be observed after waiting enough time 
for a significant flavour change to have occurred. The flavour change means changing 
probability amplitudes for the mass cigcnstatcs which corresponds to changing internal 
states. In principle, the full internal state can be measured using quantum state 
tomography |31| after tracing out the phonon modes. However, the different mass 
eigenstates can be made to have different fluorescence rates (with extra single qubit 
rotations) when coupled to an auxiliary level via an external laser field and therefore be 
distinguished by looking at the fluorescence level [12] . Then, the measured fluorescence 
level would oscillate in accordance with the neutrino flavour components. Exactly 
how it oscillates would depend on a particular implementation used, but it can be 
calculated theoretically and then compared to experimental results. Figure [3] shows, 
as an example, flavour oscillations of an electron neutrino created at t = 0, calculated 
from Eq. (j4|) with the tribimaximal mixing matrix. The kinetic and rest mass energies, 
shown in the caption, are chosen to correspond to experimentally viable numbers while 
obeying the ultrarelativistic condition. The oscillations should clearly be visible in real 
experiments with decoherence, as the decoherence time can be of the order of 10ms 
(see e.g. [21]). The oscillation period decreases if one moves closer to the normal 
realtivistic regime where the kinetic and the mass energies are similar, easing the 
requirement for a required decoherence time. 

5. Conclusion 

We have proposed an experimentally feasible scheme to simulate two or three 
generation neutrino oscillations using trapped ions. In 1+1 dimensions, our proposal 
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only utilizes experimentally proven techniques, allowing a controlled experimental 
observation of neutrino oscillations. In this work, due to its relevance to neutrino 
oscillations experiments, we have assumed that an initial electron neutrino is in an 
energy eigenstate with a momentum distribution. However, other initial conditions can 
also produce neutrino oscillations and could provide interesting alternative scenarios, 
e.g. outside the ultrarelativistic regime or initial states that are not energy eigenstates. 
Also, the current setup allows simulations of different types of neutrino oscillations 
experiments in a single experimental setup. 
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